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SUPPRESSED STAR FORMATION IN CIRCUMNUCLEAR REGIONS IN SEYFERT GALAXIES 
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ABSTRACT 

Feedback from black hole activity is widely believed to play a key role in regulating star formation and black 
hole growth. A long-standing issue is the relation between the star formation and fueling the supermassive black 
holes in active galactic nuclei (AGNs). We compile a sample of 57 Seyfert galaxies to tackle this issue. We 
estimate the surface densities of gas and star formation rates in circumnuclear regions (CNRs). Comparing with 
the well-known Kennicutt-Schmidt (K-S) law, we find that the star formation rates in CNRs of most Seyfert 
galaxies are suppressed in this sample. Feedback is suggested to explain the suppressed star formation rates. 

Subject headings: galaxies: active — galaxies: Seyfert — galaxies: feedback 



1. INTRODUCTION 

The implications of the well-known relations between black 
hole masses and bulge magnitudes (Magorrian et al. 1998), and 
the velocity dispersions (Gebhardt et al. 2000; Frarreasse & 
Merrit 2000) show a coevolution of the black holes and their 
host galaxies. However, how do black holes know the evolu- 
tion stage of the galaxies and how to control the growth of the 
black holes are currently understood via the the feedback from 
the black hole (Silk & Rees 1998; Croton et al. 2006; Schaw- 
inski et al. 2006). Numerical simulations show two roles of 
feedback from the black hole activity: (1) modulating the star 
formation rates; (2) heating the medium and finally quenching 
the black hole activity (Di Matteo et al. 2005). The direct evi- 
dence for the presence of the feedback from active black holes 
has to be shown from observations, yet. 

The main goal of the present paper is to show one piece of 
evidence for the feedback role in active galaxies. We show 
the AGN feedback domain, where starburst should be sup- 
pressed. We find the star formation rates in Seyfert galaxies 
is significantly lower than the rates predicted by the Kennicut- 
Schmidt's law. We use the cosmological parameters Hq = 
75km s" 1 Mpc" 1 , f^M = 0.3 and S1 A = 0.7 throughout calcula- 
tions. 

2. AGN FEEDBACK DOMAIN 

When the CNR medium is optically thick, namely, the opti- 
cal depth t = K a bs£ga S > 1, where K a t, s is opacity and £ gas the 
gas surface density, the radiation from the black hole activity 
will continuously heat the medium and blow the gas away so as 
to lower the star formation rates. The condition of r = 1 yields 
a critical density 
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v Kabs/5) ' Mq P C Z , (1) 

Kabs has a mean value of 5 for the CNR medium (Semenov et 
al. 2003). This is feedback driven by AGN radiation. We note 
outflows from Seyfert active nucleus have much low kinetic 
luminosities, typically ~ 10~ (3 ~ 6) Lbo1 based on X-ray warmer 
absorbers (Blustin et al. 2005), where Lb i is the bolometric 
luminosity. Feedback from outflows could be thus neglected. 
When Egas > S^ as , the AGN radiation-driven feedback will sup- 
press the star formation. On the other hand, AGN feedback 
reaches its maximum when an AGN radiates at the Eddington 
limit Lagn = L Ed d = 1.3 x 10 38 (M./M Q )erg/s. In the case of 
L E dd < ^sfr' AGN have inefficient feedback to star formation. 



With the help of SFR = 4.5 (L, R /10 44 erg s -1 ) A^yT 1 , is 
given by using the K-S law Esfr =A£g as (Kennicutt 1998a), 

m s = 8.2 x lOX^o 4 M pc" 2 , (2) 

where Esfr = SFR/irR 2 is the surface density of the star forma- 
tion rate, A = 2.5 x 10" 4 , 7 = 1 .4, M 9 = M./10 9 M Q is the black 
hole mass and R200 =R/ 200pc the size of the circumnuclear star 
forming region. When £ gas > £^ s , the gas is so dense that the 
luminosity from star formation dominates over the AGN. We 
call 

^gL < ^]gas < S, g a S , (3) 

the AGN feedback domain as shown in Fig. 1, in which the K-S 
law is broken. 

The strong radiation pressure from the black hole accretion 
disk at Eddington limit is Pagn ~ L0 x 10~ 7 Mg/?2oodyn cm -2 , 
where Mg = M./1O 8 M . The pressure from supernovae explo- 
sion is Psn ~ eSsfrc = 2.0 x 10~ 8 SsFR.2e-3clyn cm" 2 , where 
Ssfr,2 = SsFR/lO^Qyr'kpc" 2 and e_ 3 = e/10" 3 is the ef- 
ficiency converting the mass into radiation (Thompson et al 
2005). We find P AG n > 5P SN within CNRs of radius ~ 200pc 
for typical values of the parameters of e, M, and Ssfr- This 
indicates that the radiation from AGN dominates over the local 
feedback from supernovae explosion. After an AGN switches 
on, the star formation is suppressed and then feedback from 
supernovae is further weakened. The timescale of the AGN 
feedback to the starburst regions can be estimated by fpB ~ 
Egas/ /fbCLagni where Lagn is AGN luminosity, C = AH/4tt 
is the covering factor, the thermal energy £g as f» kTM glis /m p , 
k is the Boltzmann constant, m p is the proton mass, T is the 
gas temperature, M slis = ttR 2 ^^ is the gas mass and /fb is the 
feedback efficiency. We have 

(4) 
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where E gas , 4 = S gas /1O 4 M pc" 2 , T 3 = T /2 x 10 3 K, / FB - 2 = 
/fb/10" 2 , L43 = L AGN /10 43 erg s" 1 and C . 5 = C/0.5 is the cov- 
ering factor of the CNRs. Such a short timescale indicates that 
the AGN feedback is very efficient. This is supported by a large 
fraction of the post-starburst AGNs in a large Sloan Digital Sky 
Survey sample (Kauffmann et al. 2003). The physics behind 
K-S law is not sufficiently understood (Thompson et al. 2005; 
Krumholz & McKee 2005). It is beyond the scope of the present 
paper to give a quantitative description of the suppressed star 
formation rates. Comparing Seyfert galaxies with the K-S law, 
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Table 1 The Seyfert Galaxy Sample 



Seyfert 1 


Object 


Redshift 


FWHM 


logAL A 


Ref. 


logM. 


M. 


log£ gas 


SpAH 


R 


log 


SsFR 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


3C120 


0.033 




44.17 


2 


7.74" 


0.24 


4.18 


76 


0.48 


0.85 


1.32 


IC4329 


0.016 




43.32 


2 


6.99" 


0.03 


3.70 


220 


0.24 


1.31 


1.35 


MCG-2-33-34 


0.014 


1565 


42.61 


22, 5 


6.11 


0.01 


3.13 


54 


0.21 


0.70 


1.21 


MCG-5-13-17 


0.013 


4000 


43.44 


20, 1 


7.50 


0.04 


3.92 


67 


0.19 


0.79 


1.24 



NOTE.-Table 1 is published in its entirety in the electronic edition. A portion is shown here for guidance regarding its form 
and content. 



"the blackhole mass are directly taken from Peterson et al. (2004). 
*refers to [O III] FWHM. 

c based on M. -M bu i ge relation, F01475-0740: M bu i ge = -18.80; NGC 3660: M bulge = -18.38. 

Note.-(l) source name; (2) redshift; (3) FWHM of H/3 for Seyfert Is or stellar velocity dispersion a for Seyfert 2s (in km s" 1 ); 
(4) luminosity of 5100A deduced from extrapolation of F v oc v~~°- 5 or [O III]A5()07A (in erg s~'); (5) references for columns 
(3) and (4) are given below, respectively; (6) black hole mass (inA/0); (7) accretion rate (inM0yr~'); (8) gas surface density 
(in Mqpc~ 2 ); (9) surface brightness of the 3.3 /im PAH emission feature (in unit of 10 3 'ergs s~'kpc -2 ); (10) the scale of the 
starburst regions (in kpc); (11) and (12) are the lower (Eg FR ) and upper (Eg FR ) limits of surface density of star formation 
rates, respectively (in Moyr~'kpc~ 2 ). 

Reference.-(l) NED; (2) Peterson et al. (2004); (3) Spinogilio et al. (1995); (4) Doroshenko & Terebezh (1979); (5) Kinney 
et al. (1993); (6) Nelson & Whittle (1995); (7) Dahari & Robertis (1988); (8) Lipari et al. (1991); (9) Corral et al. (2005); 
(10) Kirhakos & Steiner (1990); (11) Visvanathan & Griersmith (1977); (12) Cid Fernandes et al. (2004); (13) Gu & Huang 
(2002); (14) Kailey & Lebofsky (1988); (15) Heraudeau & Simien (1998); (16) Bassani (1999); (17) Whittle et al. (1988); (18) 
Whittle (1992); (19) Garcia-Rissmann et al. (2005); (20) Crenshaw et al. (2003); (21) Marzini et al. (2003); (22) Veron-Cetty 
et al. (2001); (23) Postman & Lauer (1995). 



a universal rule of cosmic star formation, we may get the un- 
dergoing physics in the CNRs. 

We have to point out here that the short feedback time does 
NOT mean the same timescale of the starburst. The present 
fFB means the starburst rates will be suppressed once AGN is 
triggered and make it possible for AGN and starburst coexist. 

3. APPEARANCE OF FEEDBACK IN SEYFERT GALAXIES 

For the goal to test the above scenario, we compile 57 Seyfert 
galaxies (Imanishi 2002; Imanishi 2003; Imanishi & Wada 
2004). The star formation rates in CNRs of Seyfert galaxies 
can be traced by several indicators, particularly, PAH features 
at 3.3, 6.2, 7.7, 8.6 and 1 1 .2^m, which radiate from vibration of 
PAH grains containing about 50 carbon atoms. Among the fea- 
tures, 3.3/im emission is intrinsically strong and less affected 
by broad silicate dust absorption (Imanishi 2002). We choose 
3.3^m emission as an indicator of the star formation rate. We 
convert the PAH emission into IR luminosity via Li R = 10 3 Lpah 
relation with a scatter by a factor of 2-3 for pure star formation 
(Imanishi 2002). Since some PAH grains would be destroyed 
by EUV and X-ray photons from the central engine, we have 
the lower limit of the surface density of the star formation rates 

X^ FR = 35.8Lpah,4i^200 (Mgyr-'kpc- 2 ), (5) 
by using the relation of the star formation rate and the in- 
frared luminosity (eq. 7) (Kennicutt 1998a), where Lpah.41 = 
ipAH/10 41 erg s" 1 . On the other hand, the infrared emission 
from Seyfert galaxies covers the contribution from starburst and 
reprocessing radiation from AGNs, we have the upper limit of 
the surface density of the star formation rates 

E&r = 35.8L FIR , 44 ^200 (Moyr^kpc" 2 ), (6) 
where Lfir.44 = ^FiR/10 44 erg s -1 is the observed far-IR luminos- 
ity. We take the geometric average Esfr = ( ^sfr^sfr ) an d 
the error bars correspond to Sg FR and £g FR . We have to stress 



this average only represents the central value of logarithm of 
Eg FR and Sg FR and the upper and lower limits of Esfr are the 
most important. Table 1 gives the sample of Seyfert galaxies, 
which have been observed through IRTF SpeX or Subaru IRCS 
with spatial resolution of 0.9"- 1.6". 

For Seyfert 1 galaxies, we estimate LboI = 9L5100, where 
L5100 is the luminosity at 5100A and then the accretion rate 
M. = Lboi/^c 2 , where r\ = 0.1 is the accretion efficiency. The 
black hole masses are estimated from the empirical relation of 
reverberation mapping (Kaspi et al. 2000), or directly taken 
from the mapping (Peterson et al. 2004). We assume that the 
potential of the total mass within the CNRs controls the massive 
disk fueling to the black hole, where star formations are taking 
place either. Assuming the Keplerian rotation, the surface den- 
sity of the disk is 

S tot = 2. 1 x 10 V/' ^iV: 1/ X'X)o 5 M Q pc" 2 , (7) 
given by the disk model (King et al. 2002; Yi & Black- 
man 1994; Tan 2005), where the opacity K- d \, s = 5 in this re- 
gion, /. is the ratio of the black hole mass to the total, M.j = 
M./1.0M Q yr _1 and ao.i = a/0.1 is the viscosity (Shakura & 
Sunyaev 1973). This estimation is the lower limit since we re- 
place the infalling mass rates in CNRs by black hole accretion 
rates. The gas surface density of the disk £ gas = / g £ t ot 

£ gas = 1 .0 x 10 5 /g,o.05ao 4 i /5 ^i 5 A 1/5 Ml /5 R^ 5 M q pc" 2 , 

(8) 

where / g ,o.o5 = / g /0.05 is the gas fraction to the to- 
tal. Considering the disk is located inside the 
bulge, we have / g = M gas /M disk > M gas /M Bu i g e = 
(M gas /M dust ) (M dust /M.) (M. /MBuige), where M disk is the to- 
tal mass of the disk, M gas /M dust is the gas-to-dust mass ra- 
tio and MBui ge ~ 10 3 M. is the bulge mass (Kormendy & 
Gebhardt 2001; McLure & Dunlop 2002). It has been 
found that the dust mass in PG quasars is comparable 
with in Seyfert galaxies (Spinoglio et al. 2002; Haas et 
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1 2 3 4 5 6 

lo § ( S gas/^oP c " 2 ) 

FIG. 1 . — The plot of gas and star formation rate surface densities. The yellow region is the AGN feedback domain given by E gas < Eg as < Sg| s - The Compton 
thick region has Sgas > 8.0 X 1O 3 M0 pc~ 2 (i.e. Nn > 10 24 cnT 2 ). The red squares are starburst galaxies taken from Kennicutt (1998b). The cyan and blue-magenta 
stars are Seyfert 1 and 2 galaxies, respectively. The blue star is NGC 3227, in which the star formation rate is 0.05Mq yr~' and the gas mass M gas = (2 - 20) X 10 8 Mq 
within 65pc taken from Davies et al. (2006). 



al. 2003). The mean value of gas-to-dust mass ratio is 
(M gas /M dust ) ~ 250 (Haas et al. 2003). We estimated dust 
mass from M dust ~ 4.78/100^ [exp(143.38/T dust )- 1] M , 
and the dust temperature is estimated by 7d ust = (1 + 
z) [0.5-82/ln(0.3/6o M //ioo/i)] K . where f mit and / 60jU are the 
fluxes at 100/im and 60/xm in unit of Jy, respectively, D L is the 
luminosity distance in unit of Mpc (Evans et al. 2005). We find 
the mean value of (M dust /M.) = 0.2 ±0.2 in our sample. So 
we have / g > 0.05 as a lower limit in this paper. Thompson 

et al. (2005) used / g = 0.1. We note £ gas oc /. l ^ 5 , resulting in 
uncertainties of £ gas by a factor of 4 for /. = 10~ 3 - 1. a = 0.1 
is used for all Seyfert galaxies. 

For Seyfert 2 galaxies, dusty tori obscure the active regions. 
We estimate the bolometric luminosity from LboI = 3500L[o mj 
with a mean uncertainty of 0.38 dex (Heckman et al. 2004), 
where L[o mj is the [O lll]A5007 luminosity, and hence the 
black hole accretion rates. The black hole masses are estimated 
through the M, - a relation (Tremaine et al. 2002), where the 
dispersion velocity a = FWHM([0 lll])/2.35 if the dispersion 
velocity is not available. 

Fig. 1 shows the £ gas - £ S fr plot of Seyfert CNRs. We 
find that CNR gas surface densities of Seyfert galaxies are lo- 
cated within the AGN feedback domain. There are clearly three 
branches in the figure, separating the Seyfert galaxies, when 
> Sg^ > £ gas - Seyfert galaxies marked in Zone I still sat- 
isfy the K-S law. Those (Mrk 273, Mrk 938, NGC 5135 and 
NGC 1068) marked in Zone II are located between the K-S 
law and Zone III. These are ultra-luminous infrared galaxies, 
or mixed with strong starbursts. The main energy sources in 
the CNRs are in a transition state from a starburst to an AGN 
in these galaxies. The fraction of the transiting galaxies is only 
4/57 ~ 1/10. Though the completeness of the present sample is 
uncertain, this fraction implies that the transition is quite short 



and indicated by the feedback timescale from equation (4). The 
Seyfert galaxies in Zone III are undergoing suppressed star for- 
mation strongly, being 1 -2 orders lower than that predicted by 
the K-S law. The suppressed Ssfr is obviously caused by the 
feedback. Galaxies obeying the K-S law are powered by nu- 
clear energy from stars, however gravitational energy released 
from accretion onto the black holes is powering AGNs if a tran- 
sition from starburts to active galaxies happens. With the dis- 
sipation of CNR gas due to star formation and accretion onto 
the black holes, S gas is decreasing and the galaxies may return 
to the K-S law once AGNs switch off. Such a behavior likes 
evolution of stellar energy sources in the Hertzprung-Russell 
diagram. 

It has been found that black hole duty cycles follow the his- 
tory of star formation rate density (Wang et al. 2006). The 
above scenario then implies that both the black hole activities 
and starbursts are episodic (Davies et al. 2006). The multiple 
cycles of the black holes and starbursts make it impossible to 
measure the time delay between the two episodes. However 
the stellar synthesis may tell the star formation history and then 
give the black hole activity history. 

4. CONCLUSIONS AND DISCUSSIONS 

We find direct evidence for the feedback from active black 
holes in Seyfert galaxies. Once a black hole is triggered, the 
feedback will significantly suppress the starbursts within a quite 
short timescale of a few 10 5 years. The duty cycles of Seyfert 
galaxies strongly indicate there is an efficient way to frequently 
trigger black holes and quench starbursts. 

The data presented in this paper are only lower limits of the 
gas densities. Future VLT (Very Large Telescope) and ALMA 
(Atacama Large Millimiter Array) measurements of the star 
formation rates and gas densities will finally identify roles of 
the feedback from the black hole activities. 
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Seyfert 1 



Object 


Redshift 


FWHM 


log \L\ 


Ref. 


logM. 


M. 


logS gas 


5pAH 


R 


log 


^SFR 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


3C120 


0.033 




44.17 


2 


7.74" 


0.24 


4.18 


76 


0.48 


0.85 


1.32 


IC4329 


0.016 




43.32 


2 


6.99" 


0.03 


3.70 


220 


0.24 


1.31 


1.35 


MCG-2-33-34 


0.014 


1565 


42.61 


22, 5 


6.11 


0.01 


3.13 


54 


0.21 


0.70 


1.21 


MCG-5-13-17 


0.013 


4000 


43.44 


20, 1 


7.50 


0.04 


3.92 


67 


0.19 


0.79 


1.24 


Mrk79 


0.022 




43.72 


2 


7.72" 


0.08 


4.00 


71 


0.32 


0.82 


1.31 


Mrk335 


0.025 




43.86 


2 


7.15" 


0.14 


3.94 


36 


0.37 


0.52 


0.68 


Mrk509 


0.035 




44.28 


2 


8.15" 


0.30 


4.31 


62 


0.51 


0.76 


1.24 


Mrk53() 


0.029 


6560 


44.04 


20, 1 


8.35 


0.17 


4.25 


67 


0.42 


0.79 


1.15 


Mrk618 


0.035 


3018 


44.00 


21,4 


7.65 


0.16 


4.04 


44 


0.51 


0.61 


1.59 


Mrk704 


0.030 


5684 


43.53 


21, 1 


7.88 


0.05 


3.84 


32 


0.44 


0.47 


0.76 


Mrk817 


0.031 




43.82 


2 


7.69° 


0.11 


3.97 


78 


0.45 


0.86 


1.42 


Mrkl239 


0.019 


1075 


43.84 


22, 3 


6.65 


0.11 


4.02 


47 


0.17 


0.64 


1.70 


NGC863 


0.027 




43.81 


2 


7.68° 


0.10 


4.00 


33 


0.40 


0.48 


0.95 


NGC931 


0.016 


1830 
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"the blackhole mass are directly taken from Peterson et al. (2004). 
^refers to [0 III] FWHM. 

c based on M. -M blllge relation, F01475-0740: M bulge = -18.80; NGC 3660: M bulgc = -18.38. 

Note.-(l) source name; (2) redshift; (3) FWHM of H/3 for Seyfert Is or stellar velocity dispersion a for Seyfert 2s (in km s~'); 
(4) luminosity of 5100A deduced from extrapolation of F„ oc v~° 5 or [O III]A5007A (in erg s~'); (5) references for columns 
(3) and (4) are given below, respectively; (6) black hole mass (in Mq); (7) accretion rate (in M^y^T 1 ); (8) gas surface density 
(in Mqpc~ 2 ); (9) surface brightness of the 3.3 fim PAH emission feature (in X 10 39 ergs s~'kpc~ 2 ); (10) the scale of the starburst 
regions (in kpc); (11) and (12) are the lower (Sj FR ) and upper (Sj FR ) limits of surface density of star formation rates, respectively 
(in M0yr~'kpc~ 2 ). 

Reference.-(l) NED; (2) Peterson et al. (2004); (3) Spinogilio et al. (1995); (4) Doroshenko & Terebezh (1979); (5) Kinney et 
al. (1993); (6) Nelson & Whittle (1995); (7) Dahari & Robertis (1988); (8) Lipari et al. (1991); (9) Corral et al. (2005); (10) 
Kirhakos & Steiner (1990); (11) Visvanathan & Griersmith (1977); (12) Cid Femandes et al. (2004); (13) Gu & Huang (2002); 
(14) Kailey & Lebofsky (1988); (15) Heraudeau & Simien (1998); (16) Bassani (1999); (17) Whittle et al. (1988); (18) Whittle 
(1992); (19) Garcia-Rissmann et al. (2005); (20) Crenshaw et al. (2003); (21) Marzini et al. (2003); (22) Veron-Cetty et al. 
(2001); (23) Postman & Lauer (1995). 
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